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We provide a comprehensive picture of magnetotransport in graphene monolayers in the limit of 
non-quantizing magnetic fields. We discuss the effects of two carrier transport, weak localization, 
weak anti-localization, and strong localization for graphene devices of various mobilities, through 
theory, experiments and numerical simulations. In particular, we observe the weak localization of the 
localization length, which allows us to make the connection between weak and strong localization. 
It provides a unified framework for both localizations, which explains the observed experimental 
features. We compare these results to numerical simulation and find a remarkable agreement between 
theory, experiment and numerics. Various graphene devices were used in this study, including 
graphene on different substrates, such as glass and silicon, as well as low and high mobility devices. 



I. INTRODUCTION 

Graphene has attracted a considerable amount of at- 
tention due to the ease in isolating a single sheet of 
graphite via mechanical exfoliation.^^ Despite the fact 
that it is only one atom thick, exfoliated graphene has 
shown extraordinary transport properties and can be 
used as a novel material for many potential applications. 
Its unique band structure has led to many interesting^ 
phenomena such as tunable charge c arri ers densitieSj- 
anomalous integer quantum Hall efFeclP^ and ultrahigh 
mobilities at room temperature.^ The most noteworthy 
property of the band structure is the existence of two 
degenerate Dirac cones,!^ which leads to two degenerate 
valleys (K and K'). The existence of these valleys, each 
with a linear dispersion, is the main reason that graphene 
has transport properties which differ from most other 
semiconductors or semimetals. At high magnetic fields, 
a striking example is the quantum Hall effect, where the 
Hall conductance quantization occurs in steps of Ae^/h, 
because of the spin and valley degeneracy, and where the 
Landau level quantization is proportional to V~B, where 
B is the perpendicular magnetic field component. This 
square root dependence leads to a very large lowest Lan- 
dau level splitting, where the field quantization has been 
observed up to room temperature.13 While the four-fold 
degeneracy can be lifted at very high fields in high mo- 
bility graphene, it is invisible at lower fields, which is the 
focus of this work. 

The goal of this work is to present a comprehensive 
understanding of the magnetoresistance of monolayer 
graphene in a regime where the Landau quantization is 
not important, i.e., fiqB <^ 1, where Hq is the quan- 
tum mobility. In this regime, the Landau quantization 
does not play any role, but the valley degeneracy, the lin- 
ear dispersion and the underlying disorder potential lead 
to interesting magnetotransport phenomena such as two 
carrier transport (2CT), weak localization (WL), weak 
anti-localization (WAL), and strong localization (SL). 

Most of the magnetotransport properties of graphene 
can be understood in terms of two types of scattering 
mechanisms: intervalley scattering, where electrons are 
scattered from one valley to the other (K to K') with 



a rate hereby noted as t^^ and intravalley scattering, 
where electrons scatter within a valley as described by a 
rate t~^. In general, intervalley stems from short range 
scattering, such as atomic defects, including grain bound- 
aries, whereas intravalley scattering is long range and is 
typically stronger and also includes large scale inhomo- 
geneities, as well as charged impurities in the substrate. 
Long range potential variations are responsible for the ex- 
istence of both carriers types (electrons and holes) simul- 
taneously at the charge neutrality point (CNP). The car- 
rier distribution can be understood by assuming that the 
total carrier density is given by n+p — ^/n-Q + {n — p)"^^ 
where n—p — C|Vg— Vol/e, C is the gate capacitance, Vg 
the gate voltage, and Vd the gate voltage at the charge 
neutrality point (CNP) with a total residual density no- 
n and p are the carrier densities of electrons and holes, 
respectively. This allows us to define electron and hole 
mobilities, /j,„ = cr„/en and /ip = cTp/ep, where ct„ and ap 
are the corresponding electron and hole conductivities. 
Experimentally, /i„ and /ip are often slightly different, 
which is likely due to the asymmetry of the scattering 
potential for electrons and holes. 



II. TWO CARRIER TRANSPORT 

We start our magnetotransport analysis, by discussing 
the simplest non-trivial contribution, which is due to the 
simultaneous existence of the two types of carriers near 
the CNP. As long as jigB <^ 1, this contribution can 
be evaluated using Drude's expression for two carriers, 
where p^^ = a^^ and p^y = sign{qri) B / rj for the two 
types of carrier densities [rj = n or p). The total resis- 
tivity is then given by ptot = {^n + o'p)^^, where are 
the single band conductivity matrices. This yields the 
following field dependence of the magnetoresistivity 
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From this it follows that the relative field dependence can 
be written as 
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where 
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Expressing the conductivities in equation ([s]) in terms 
of mobilities, equating the zero field resistivity to pxx = 
((T„ + <Jp)~^ and identifying Pmax as the zero field resis- 
tivity at the CNP, we obtain 
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Hence, as a function of B the 2CT gives rise to a 
parabolic positive relative magnetoresistance before sat- 
urating to n()/\n — p\ when i? ~ _Bi, which is valid as 
long as HqB ^ 1. Close to the CNP, i.e., |n — p| ^ ng, 
Bq is simply equal to the inverse median mobility and 
Bi ^ Bq, which is the regime where the 2CT effect is 
the largest. Away from the CNP ^ uq), we have 

a maximum Ap/p ~ no/fiH, which decreases away from 
the CNP, since uh — \n — p\, which is the Hall density, 
increases. The 2CT will play an important role, as long 
as pqB <C 1, to give rise to a positive magnetoresistance 
due to large scale inhomogeneities. This effect is very 
important in graphene as compared to other two dimen- 
sional systems because of the absence of a gap between 
the electron and hole carriers. The 2CT will serve as ba- 
sis for the understanding of the classical contribution due 
to long range scattering, which will be present at all tem- 
peratures and even increases with temperature, since no 
increases with temperature due to the thermal activation 
of the electron and hole carriers. 



III. THEORY OF WEAK LOCALIZATION 



second term represents the short time cut-off (rg), below 
which no elastic scattering occurs, and the third term 
is the phase coherence time (r^) cut-off, beyond which 
phase coherence is lost. At low temperatures, where 3> 
Te, this leads to the typical logarithmic WL correction of 
the conductivity Sa = ~^^Pret ^ ln{T^/Te). 

To evaluate the magnetic field dependence of the re- 
turn probability Psit), we have to solve for the return 
probability according to the field dependent diffusion 



equation^ 



+ D{i^ + '^-^Af]PB{r, /, t) - 5{r - r')5{t), (6) 

where A is the vector potential. The solution to this 
diffusion equation can be evaluated for r — r' and is 
given by 



PB{t) = 
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where tb = Ti/AeBD. At zero field we recover 
P[){t). In the long time limit, where we have Psit) — 
(ei?/ft-)e^*/^'^^ , this leads to the field induced destruc- 
tion of the phase coherence when approximately one flux 
penetrates the area = Dtb- The full magnetic field 
dependence is obtained by inserting PB{t) into the ex- 
pression of the coherent return probability and evaluat- 
ing the integral, i.e.. 
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where B^ — h/AeDr^ and B^ = Ti/AeDTf. are the char- 
acteristic magnetic fields for and Te, whereas ip is the 
digamma function. The magnetic field correction to the 
coherent return probability, is then given by 



Moving beyond the Drude and classical description of 
transport, we need to include quantum effects, of which 
coherent backscattering is the most important contribu- 
tion. Coherent backscattering leads to WL .121121 To obtain 
an expression for the WL correction, we have to evaluate 
the return probability of all possible trajectories'^^ At 
zero magnetic field, the coherent return probability P^et 
can be expressed asl^^ 
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where the first term {AnDt) ^ in the integrant is the 
return probability {Po{t)) for diffusion constant D, the 
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where F{z) = ip{l/2 + 1/z) + ln(z), F{z) ~ z^/24: 
for z < 1, F{z) ~ ln(l + z/4e'^) for z > 1, and 7 
is the Euler constant. The second term in equation 
([9]) is responsible for the reduction of the return prob- 
ability due to the presence of a magnetic field, which 
leads to the observed peak in resistance of width ~ B^ 
around zero magnetic field, characteristic of weak local- 
ization in disordered systems.l^ The first term renormal- 
izes the coherent return probability at large fields, i.e., 
inDAPretiB ^ 00) ^ \n{B^/{B^ + Be)). This ex- 
pression leads to the usual WL correction to the con- 
ductivity when expressing the conductivity correction in 
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terms of the coherent retur n probabihty correction, i.e., 
Act = -Ae'DAPret/h^^ 

In deriving equation (|9| we assumed the existence of 
only two types of scattering mechanisms (phase break- 
ing) and Tg (elastic). However, in graphene, we have 
two very different types of elastic scattering mechanisms: 
those that scatter electrons within a valley (intravalley. 
To) and those that scatter between valleys (intervalley, K 
to K', Te). For coherent backscattering, only intervalley 
scattering contributes, which also includes inelastic scat- 
tering (r^). Hence, the short time cutoff scattering time 
is provided by Te in graphene, which is short ranged and 
can be due to grain boundaries or lattice defects. How- 
ever, intravalley scattering is usually much stronger in 
graphene, but doesn't contribute to coherent backscat- 
tering and can lead to enhanced forward scattering at 
zero magnetic field and positive magnetoresistance. 

In an important work by McCann and co-workers,^ 
based on earlier work in honeycomb lattices,!^ the au- 
thors have obtained a general expression for the WL and 
WAL correction specific to graphene, which determines 
the dependence of the magnetoresistivity as a function of 
B involving the scattering parameters Tq and Te explic- 
itly. They obtained the following expression: 
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FIG. 1: The relative return probability is shown as a function 
of the magnetic field for the different models. The green curve 
is from expression ^ , the blue curve is from expression 1 10 1 
with = SO-B^, and the red curve is obtained from expres- 
sion ( |19[ l with Bo = 5005^ and B\ = co. For all curves we 
assumed B, = Bs- 
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where Bi, = Ba + Bf,/2. This expression can also be ob- 
tained directly using equation (|9| by replacing the short 
time cut-off 6"*/"^= in the magnetic field correction by 
(■g-t/r^ -f 26"*/"^*). While the first two terms in equation 



( 10 1 lead to the usual WL localization effect as in equa- 
tion ([9|, the last term leads to WAL due to the pres- 
ence of intravalley scattering as illustrated in figure 1. 
This term becomes important when B approaches Bj, 
and leads to a positive magnetoresistance contribution. 



IV. STRONG LOCALIZATION 

In some graphene samples disorder can be very im- 
portant, like for instance in intentionally disordered 
exfoliated graphene, which was shown to lead to 
strong localization.'^^ Strong localization or Anderson 
localizatiorfi^ is obtained when the transmission is ex- 
ponentially suppressed, i.e., when a coherent wave is 
backscattered due to wave interference. Hence, strong 
localization is a direct consequence of coherent backscat- 
tering. In the language of Anderson localization, the im- 
portant parameter is the localization length Lc^ which 
measures the exponential increase of the resistance with 
size. In a system where the coherence length is infinite, 
the inverse localization length can be defined by 
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where L is the length of the system and R the resistance. 
In order to incorporate strong localization in our discus- 
sion on weak localization, we can define the localization 
time as Tc — L^/D, which represents the time before 
the charge carrier is localized in a diffusive system. The 
coherent return probability then has to be modified ac- 
cording to 
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where Be = h/AeDTc = H/AeLl and Bc^, ^ Bc + B^. The 
effect of strong localization, is therefore simply to replace 
B^ by -Bc0- This is very intuitive, since if Lc L^, Lc 
will play the role of in determining the magnetic field 
correction, since no closed trajectory can occur beyond 
Lc- The relative change in coherent return probability 
can then be written as 
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where the field dependence is mainly determined by B^^. 
Since we define the localization length assuming a co- 
herent system, the field induced relative change in Lc is 
directly determined by the change in the relative coher- 
ent return probability for zero dephasing: 
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VI. NUMERICAL SIMULATIONS 
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Equation ( 14 ) tells us how strong localization is con- 



nected to the coherent return probability and therefore to 
WL. Since it is quite challenging to obtain an expression 
for the localization length in two dimensions directly, we 
will resort to evaluating the localization length numeri- 
cally below. 



V. GRAPHENE NANORIBBONS 

In many situations, both experimentally and numeri- 
cally, one deals with a situation where the sample width is 
finite. This leads to a different field dependence, since the 
sample effectively becomes quasi-one dimensional (QID) 
instead of 2D, when > W and W is the width of the 
graphene ribbon. In this limit, the governing diffusion 
equation is no longer 2D, but has to be replaced by the 
ID analogue, i.e., Psit) = (l/V47rDi)e-*/^\E^ where the 
first term represents ID diffusion at zero field and the sec- 
ond term represents the field induced destruction of phase 
coherence in QID. is given as usual by L^/D, where 
BWLb = y/sjih/e scales with the flux quantum through 
the relevant characteristic area WLb^ This leads to 



(l_e-*/^=)e-*/^*dt 
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When ^ + Tf^^ we obtain the following ex- 
pression for the magnetic field correction of the coherent 
return probability: 
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Recalling that Act = -4e'^DAPret/h, we recover the 
expression for the WL correction to the conductivity ob- 
tained by McCann and co- workers. 

For the strong disorder case, we have to again intro- 
duce the localization time Tc, which has the effect of sim- 
ply replacing l/L^ by l/Ll^ = l/i| + l/Ll. We then 
obtain for the relative change in the coherent return prob- 
ability. 
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where Bw = fi/^eW'^. This is the result for the QID 
limit, when or 3> W. Here again we have 
ALc/Lc — —APretl Pret\B^=n relative change of 

the localization length. 



We now move to evaluate numerically for graphene 
assuming ^ L^j W . There are two different limits: 
2D and QID, i.e., Lc (the high disorder limit) and 

W (the low disorder limit), respectively. These two 
limits are mainly determined by the amount of disorder 
in the system. For the numerical simulations, we will only 
consider short ranged disorder, since this is the relevant 
source for WL. The effect of long range disorder is to 
induce the co-existence of electrons and holes close to 
the charge neutrality point as well as WAL. 

To model the system accurately, we consider a tight 
binding model with a honeycomb structure with two 
types of edges: armchair and zigzag. The hopping term 
is given by t ~3 eV. The two ends of the graphene device 
are assumed connected to wide ohmic contacts with a 
quadratic dispersion. The short range disorder is imple- 
mented by adding a random onsite potential Vi on every 
atom site uniformly distributed with —V/2 < Vi < V/2. 
The magnetic field is simply a phase factor (j) in the 
hopping term, where 27r0 = 1 corresponds to one flux 
quantum in a hexagon. The two terminal transmission 
probability T is then evaluated iteratively as a function 
of the length, L, of the system, using an efficient itera- 
tive Green's function method.^^ The average two termi- 
nal conductance {T{L))e'^ /h, resistance {l/T{L))h/e^ or 
logarithm conductance (ln(T(L))) are obtained by aver- 
aging over many disorder configurations ((•)). For a given 
W and at large enough L, the system is always localized 
due to Anderson localization and Lc can be extracted 
whenever L ^ Lc- 

We now move to compare the numerically obtained Lc 
with the analytical expressions derived above for the rel- 
ative change in Lc, i.e., testing equation (14). In the 2D 



limit, where W ^ Lc, the field dependence of the coher- 



ent return probability is determined by equation (13 1, 
where Bc^ — Be in the coherent limit. The numerical 
field dependence of the relative Lc is then rescaled by 
B /Be, where Be ^ L~^ and can be compared to the rel- 
ative coherent return probabilities as shown in figure [2] 
The agreement is quite remarkable with only two fitting 
parameters Be and Be- Away from the Dirac point we 
expect to have BejBc — W, where W is approximately 
the number of quantum channels in a quasi one dimen- 
sional system. We indeed obtained Be/ Be — 100 for the 
comparison to the numerical data where 40 < Vl^ < 120 
in units of the lattice constant. 

In the QID limit, where W <C Lc, the field dependence 
is rescaled by B/^/l2BwBe, where \/BwBc 



determined from equation (17) 



l/WLe 



In both limits, the agreement between the numerically 
determined relative localization length and the relative 
coherent return probability is excellent, which confirms 
that indeed ALc/Lc — —APret/ F'ret\B^=o- Hence, the 
field dependence of the coherent return probability de- 
termines Be, which in turn determines Lc and therefore 
provides us with a way to extract the localization length 
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FIG. 2: We compare the a naly tical expression (solid line) for 
APret / Pret from equation ( 13 1 with the numerically obtained 



relative change in the localization length using the parameter 
Be = lOOPc- The different dot colors (green, blue, red) rep- 
resent different widths (40 x a, 80 x a, 120 x a) of the simulated 
graphene device, where a is the lattice constant. 



FIG. 3: We compare the analytical expression (solid line) for 
l^Pret/Pret from equatlou (17 1, relevant to the quasi ID sit- 



uation, with the numerically obtained relative change in the 
localization length using the parameter B^ = lOOPc- The dif- 
ferent dot colors (blue, red) represent different widths (80 x a, 
120 X a) of the simulated graphene device, where a is the lat- 
tice constant. 



simply from the magnetic field dependence. 

It is now possible to make the connection between Lc 
and the resistance, since hi{Rxx) — > L/Lc- Therefore, 
small variations in Lj, lead to AR^j./Rxx ~ —LALc/L^- 
This is assuming that ^ L. However, in experiments, 
where this is often not the case, then plays the role 
of effective sample size (L^ — L) so that when <^ L 
we can write AR^x/ Rxx — —L^AL^/L'^ instead. This is 
often used in experiments to extract Lc from the temper- 
ature dependence of the resistance, since also depends 
on temperature.!^ Here, we restrict ourselves to a fixed 
temperature, where is constant, but is magnetic 
field dependent. Hence, when <^ L wc have using 
equation (14) 
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For small variations {AR/R <^ 1), different contribu- 
tions to the resistance are additive, hence the total effect 
on the resistance, including the two carrier effect, yields 
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, Lc \n{l + Bc/Bc^) 1 + (S/Bi)2' 

(19) 

which allows the extraction of all relevant length scales 
from the magnetic field dependence alone, recalling that 
L~'^ + L^'^ = 4eBc^/h, = 4eBe/h, Bq ~ Pmax/ iJ-Pxx, 
and Bi ~ BqUq/uh- 



VII. EXPERIMENTS: 



While WL was observed in graphene in a number of 
previous experiments,'2^'^ and fitted using McCann's 
WL expression,'^] they were generally limited to low dis- 
order. Here we will use expression ( 19 ) derived above 



and valid for all disorder strengths and fit it to the exper- 
imental data. We performed experiments on large scale 
graphene as well as lithographically defined Hall bars and 
graphene nano-ribbons, in addition to large (over 100 
fim) single crystal graphene. Monolayers of graphene 
were grown by chemical vapor deposition (CVD) of hy- 
drocarbons on 25 /xm-thick commercial Cu foils. The 
CVD process used was similar to those described in pre- 
vious work^SlHlIl 

We start with low mobility samples, which cannot be 
fitted with the standard WL theorySI alone because of 
the importance of strong localization in this regime. In- 
stead we use the expression derived in equation (13 1, 
which incorporates the effect of strong localization. In 
figure |4] we show the relative resistance change of a rep- 
resentative low mobility graphene Hall bar. Typical fea- 
tures of low mobility samples are a wide peak around 
zero magnetic field and an important relative resistance 
change. Using the best fit over the entire magnetic field 
range, we obtain L^ ~ 50nm and Lc — 120nm, which 
shows that the localization length is of the same order of 
magnitude as the phase coherence length. The agreement 
between the fit and the data is quite remarkable over the 
entire available magnetic field range and only relies on 
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FIG. 4: The measured magneto-resistance of a low mobility 
(/X < lOOcm^/Vs) graphene Hall bar device at lOOmK. Shown 
is the relative resistance (blues dots) and the fit using equation 
( [T9| with fitting parameters B^^ — 80mT, = lOT and 
L^/Lc = 0.4 (red line). Bo and B\ are larger than lOOT and 
do not affect the fit. 



three fitting parameters, in this case Sc^, Be and L^/ L^. 
The other parameters -Bq and Bi are too large to signifi- 
cantly affect the fit. For even higher disorder, we expect 
Lc to be smaller than and we would need to mod- 
ify our assumption of the linear approximation used in 
equation (18). Instead, we expect the change in mag- 



netoresistance to depend exponentially on Lc{B) and its 
associated field dependence. This is amplified further for 
long but narrow samples (nanoribbons). 
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FIG. 5: The measured magneto-resistance of a large scale 
graphene sample on glass (blue dots). Shown is the relative 
resistance and the fit (red line) using equation ( |19[ ) with fitting 
parameters Bc4, = 3.5mT, L^/Lc = 0.0925, Bo = 55. IT and 
Be = 9T. 



The next step is to look at a large scale (~cm) 
graphene sample grown by CVD deposited on a glass 
slide with evaporated gold as contacts. The magnetic 
field dependence is shown in figure [5] which shows a rel- 
atively broad peak at zero magnetic field. There is also 
a parabolic increase in resistance at large field, which we 
attribute to 2CT. Using equation ( 19 ) to fit the data 
we extract the following length scales ~ 220nm, 
Lc — 2.3/xm, and ~ 4.3nm. Here again, the agree- 
ment between the fit and the data is quite remarkable, 
considering that the fit extends over the entire magnetic 
field range. 
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FIG. 6: The measured magneto-resistance of a large scale 
graphene sample on Si02/Si (blue dots). Shown is the sym- 
metrized (Rxx{B) -\- Rxx{—B)) relative resistance and the 
fit (red line) using equation 
Be0 = 1.4mT, L^/Lc = 0.055 



(191 



with fitting parameters 
Bo = 11. IT, Be = llmT 



and Bi — 20. ST and a zoom-in of the low field part in the 
inset. 

The next sample we are considering is a large scale 
(^cm) sized graphene sample grown by CVD and de- 
posited on a standard Si02/Si substrate, where the 
doped silicon can be used as a back gate. The main differ- 
ence is that this sample was grown using isotopically pure 
C^^ methane gas instead of C^^, so that the Raman peaks 
are shifted by ^/l2/13^ However, the different isotope 
does not affect the magnetotransport. The magnetic field 
dependence is shown in figure [Tj which shows a narrow 
weak localization peak at zero magnetic field. There is 
also a large parabolic increase in resistance at large field, 
which we attribute again to the 2CT. Using equation 
( 19 1 to fit the experimental data we obtain ~ 330nm, 
~ 6/im, Lf, ~ 120nm and a mean field effect mobility 
of fi ~0.18m^/Vs, which is common for large scale CVD 
graphene, and a residual density of no ~ 7.0 x lO^^cm"^. 

Interestingly, this is a prime example of a large 
parabolic-like positive magneto-resistance background, 
which cannot be fitted using WAL as per equation (10 1, 



since WAL gives a negatively curved magneto-resistance 
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at large fields. Here, we clearly demonstrate that the ori- 
gin of the large positive magnetoresistance background is 
likely due to the 2CT because of the coexistence of both 
n and p type carriers. Strikingly, the quality of the fit 
is very remarkable over the entire available experimental 
magnetic field range. 
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FIG. 7: The measured symmetrized magneto-resistance 
{Rxx{B) + Rxx{—B)) of a large scale graphene sample on glass 
measured at Vg ~ -1-801^ (blue dots) . Shown is the relative re- 
sistance and the fit (red line) using equation ( 19 1 with fitting 
parameters B^^, = 0.25mT, L^/Lc = 0.0067, Bo = 18. 7r, 
and Be = 3.5mT 

We now turn to the higher mobility case, where we con- 
sider a graphene field effect device made out of a large 
(~ 200/im diameter) single graphene crystal grown by 
CVD. Because of its single crystal nature, this sample 
has a mobility comparable to exfoliated graphene on sil- 
icon oxide. From the gate voltage dependence we ex- 
tracted a field effect mobility of about /i„ ~ 0.43m^/Vs 
for electrons and fip ~ 0.63m^/Vs for holes. The CNP 
was at Vg — 5AV. Using equation (19) to fit the ex- 



perimental data we can extract the various length scales 
and obtain ~ 810nm, Lc — 0.1mm, Le — 220nm 
and ^yjLnJip — 0.5. For the median mobility, we used 
the experimental resistivity ratio of Pmax/p{+^OV) — 9 
and find a remarkable agreement to the field effect mo- 



bility extracted from the gate voltage dependence. For 
the residual density we obtain no ~ 5.8 x lO^^cm"^ using 

fl^^ = enoPmax- 

The overall experimental behavior can be summarized 
as follows: at low temperatures, very low mobility sam- 
ples show a very wide peak in the resistance at zero field 
(as shown in figure|4|. In this regime it is important to in- 
clude strong localization effects to understand the width 
of the peak, since Lc is comparable to L^. The effective 
width is then determined by ~ l/^c + 1/-^^- With in- 
creasing mobility, the peak becomes increasingly sharper 
(figures 5-7) and the relevant length scale becomes L^, 
who's inverse square determines the peak width. In addi- 
tion to the WL effect (sharp negative magnetoresistance) , 
which all samples show, there is often a parabolic posi- 
tive magnetoresistance which is mainly due to the 2CT 
effect because of large scale inhomogeneities in addition 
to WAL in very clean samples. 



VIII. CONCLUSION: 

We presented a comprehensive picture of the observed 
magnetoresistance in graphene spanning all disorder lev- 
els. At high disorder, localization becomes the dominant 
mechanism, where the field dependence can be under- 
stood in terms of the field dependent localization length, 
which interestingly follows the field dependence of the co- 
herent return probability and hence the field dependence 
of the the WL correction. In this picture, WL is simply 
a consequence of the field dependent localization length. 
Both, WL and strong localization, are due to short-range 
scattering. Long range scattering causes a positive mag- 
netoresistance effect on top of the negative WL magne- 
toresistance, which stems from WAL (due to intravalley 
scattering) and from the presence of two types of carriers 
close to the CNP (the 2CT). While WL is generic to all 
disordered systems, both WAL and the 2CT are specific 
to graphene, where both play an important role in the 
understanding of magneto-transport. 
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